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ABSTRACT 

Using a sample of 24264 optically selected AGNs from the SDSS DR7 database, we char- 
acterise how the profile of the [O lll]A,5007 emission line relates to bolometric luminosity 
(Lagn), Eddington ratio, radio loudness, radio luminosity (Li.4ghz) and optical class (i.e., 
broad/narrow line Seyfert 1, Type 2) to determine what drives the kinematics of this kpc- 
scale line emitting gas. Firstly, we use spectral stacking to characterise how the average 
[O lll]A-5007 profile changes as function of these five variables. After accounting for the 
known correlation between Lagn and Li 4 ghz, we report that Li 4 ghz has the strongest in- 
fluence on the [O lll]A-5007 profile, with AGNs of moderate radio luminosity (Li.4ghz= 
10^^ - 10^^ W Hz"') having the broadest [O lll]A,5007 profiles. Conversely, we find only 
a modest change in the [O lll]A-5007 profile with increasing radio loudness and find no sig- 
nificant difference between the [O lll]A-5007 profiles of broad and narrow-line Seyfert Is. 
Similarly, only the very highest Eddington ratio AGNs (i.e., > 0.3) show any signs of having 
broadened [O lll]A,5007 profiles, although the small numbers of such extreme AGNs in our 
sample mean we cannot rule out that other processes (e.g., radio jets) are responsible for this 
broadening. The [O lll]A-5007 profiles of Type 1 and Type 2 AGNs show the same trends in 
terms of line width, but Type 1 AGNs display a much stronger "blue-wing", which we inter- 
pret as evidence of outflowing ionised gas. We perform multi- component fitting to the Hp, 
[O lll]U4959,5007, [N ll]U6548,6584 and Ha fines for all the AGNs in our sample to cal- 
culate the proportions of AGNs with broad [O lll]A,5007 profiles. The individual fits confirm 
the results from our stacked spectra; AGNs with Li 4 ghz> 10^^ W Hz"' are roughly 5 times 
more likely to have extremely broad [O lll]A,5007 lines (FWHMAvg> 1000 km s"') compared 
to lower Li 4 ghz AGNs and the width of the [O lll]A,5007 line peaks in moderate radio lu- 
minosity AGNs (L1.4 GHz"^ lO^'' W Hz"'). Our results are consistent with the most disturbed 
gas kinematics being induced by compact radio cores (rather than powerful radio jets), al- 
though broadened [O lll]A-5007 lines are also present, but much rarer, in low L1.4 ghz systems. 
Our catalogue of multi-component fits is freely available as an online resource for statistical 
studies of the kinematics and luminosities of the narrow and broad line AGN regions and the 
identification of potential targets for follow-up observations. 

Keywords: Keywords 



* Catalogue of multi-component fit data available at: 

http://sites.google.com/site/sdssalpaka 

t E-mail: j.r.mullaney@dur.ac.uk 



1 INTRODUCTION 

The most successful models of galaxy evolution invoke interactions 
between accreting supermassive black holes (hereafter, SMBHs) 
and their host galaxies to reproduce the most fundamental prop- 
erties of today's g alaxy population (e.g.. iDi Matteo et alj|200j : 
iBower et al.l |2006|) . Many of these models assume that a frac- 
tion of the energy released by an accreting SMBH (i.e., an ac- 
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live galactic nucleus; hereafter, AGN) drives gas and dust from 
the central regions of a galaxy, quenching star-formation (e. g., 
iDi Matteo et"al]|2005l : ISpringel et aLll2005l: iHopkins et 
However, it is not yet clear what process drives these galaxy-scale 
"outflows", with both radio jets (either compac t or extended, e.g., 
iNesvadba et al.l 2006l,'200 8l:lRosario et alj2010h and radiation pres- 



2OIOI) having been suggested as viable 



sure (e.g. [Alex ander et al 
mechanisms. 

Most observational evidence for galaxy-scale outflows has so- 



Nesvadba et al. 20061 2008 


; Alexander et al. 2010l; iRosario et al.l 


2010: Greene etal.ll2011; 


Harrison et al] 2012). These studies. 



while providing considerable insight into the scale and detailed 
kinematics of the outflows, are subject to strong selection biases 
(e.g., pre-selected for radio loudness/luminosity, bolometric lumi- 
nosity, unique emission line properties). Therefore, it has so-far 
been difficult to place these outflows in the context of the general 
AGN and galaxy populations and establish which AGN properties 
are most conducive to their production (i.e., radio jets/cores, high 
bolometric luminosity, high Eddington ratio etc.). In this study we 
aim to address this by taking the complementary approach of ex- 
ploring in coarser detail the gas kinematics of much larger num- 
bers of AGNs with physical properties spanning large sections of 
parameter space. 

The primary diagnostics that the above studies used to ini- 
tially identify kpc-scale outflows were broad, often blueshifted 
[O IIl]A,5007 emission lines in integrated (i.e., ID) rest-frame op- 
tical spectra. Since [O ni]A.5007 is produced through a forbidden 
transition it is only emitted by low-density gas. As such, any broad- 
ening or shifting of the [O III] ^5007 emission lines are the result of 
strong velocity gradients or bulk motions (i.e., disturbed kinemat- 
ics) in the so-called narrow-line region (here after, NLR) o f AGNs 
which can be extended over kpc-scales (e.g. |Poggelll989l) . as op- 
posed to the dense, sub-pc sized broad-line region. 

While broadened (beyond that expected from galaxy rota- 
tion) and asymmetric [O ni]A.5007 profiles in AGN spectra have 
been reported for ov e r three decades (earl y studies include, e.g. , 
Heckman etai]|l98ll: iFeldman et a"i]|l982l : iHeckman etZlll984 



Whittlelll985ah . there remains some uncertainty as to what phys- 



ical process cause the gas kinematics that result in these pro- 
files. Early studies of small to moderately sized samples of nearby 
AGNs (numbering < 150) reported positive correlations between 
[O IIl]X5007 profile widths and both [O ni]A,5007 and 1.4 GHz 
radio luminosities for flat-spectrum radio so urce s with L; 4 ghz< 
3 X 10^"^ W Hz~'(e.g..lHeclanan et alj[T984l. also I Wilson jfeWilUsI 



ll980l:IWhittlj|l985b[). However, since Ljn tt t I is broadly correlated 



with L; 4 r.HT(e.g..lBaum & H eckmanlll98 9'. Rawli ngs et alj 



19891 : 



2004 



Zirbe l & Baum 1995: Tadhu nter et all [1998: Wil ls et all 
de Vries et al.i i2007). these early studies were unable determine 
which of these parameters were more fundamentally linked to the 
kinematics of the narrow-line region. With the availability of large 
spectroscopic data sets (e.g., the SDSS and 2dF surveys), more re- 
cent studies have explored how the [O IIl]A,5007 profile changes 
across wide areas of parameter space. For example, using a sample 
of 1749 SDSS-selected AGNs. Greene & Ho (2005a) found that the 
narrow "core" of the [O IIllA,5007 line traced the kinematics of the 



' This "quasai-mode" feedback is distinct from the "radio-mode", in which 
the intracluster medium is prevented from cooling onto galaxies (where it 
can form stars) by powerful, extended, radio jets produced by a low accre- 
tion rate, radio-loud AGN. 



host galaxy, with an y departures corr elating with the E ddington ra- 
tio of the AGN (also leian et alj2006l : although see also lZhang et al.l 
I20T 11). However, while some recent studies have explored, in detail, 
the links between gas kinematics and the radio properties for small 
samples of AGNs (e.g., Wu 20091: iHusemann et alJlioTl) . analysis 
of much larger samples are needed to break the degeneracy between 
radio and [O IIl]X5007 luminosities to determine how gas kinemat- 
ics relate to the radio properties for the AGN population in general. 

In this study, we use the SDSS spectral database to derive the 
average [G IIl]A,5007 line profiles of a large sample (i.e., number- 
ing 24,267) of optically-selected AGNs at z < 0.4 and explore how 
these average profiles change as a function of key AGN parame- 
ters (specifically, [O IIl]X5007 luminosity, Eddington ratio, radio 
loudness, radio luminosity and spectral classification [i.e., optical 
Type 1, Type 2]). While the single aperture spectra provided by the 
SDSS provide no information regarding the extent of any kinemat- 
ically disturbed gas, the SDSS spectral archive provides the statis- 
tics needed to determine how gas kinematics relate to other AGN 
properties (including the radio properties) and establish the preva- 
lence of high gas velocities among different AGN sub-populations. 
Furthermore, it is important to stress that the redshift range of our 
sample does not span th e epoch of peak star-formation and AGN 
activi t y (i.e., z ^ ^ ~ 3 ; iHopkins & Beacon] I2OO6I : iMerloni et alj 
I2OO4I : lAird et al. I I20T0I) when quasar-mode AGN driven outflows 
are thought to have had their greatest impact on galaxy growth. 
However, the high signal-to-noise ratios and large statistics allowed 
by surveys of the local Universe can give us a better insight into 
which physical processes are most capable of disturbing the gas 
reservoirs around an AGN. 

In the following section we outline our initial source selection 
procedure, describe our multi-component line fitting routine and 
outline how the results from this routine were used to reclassify 
and analyse the AGNs in our sample. In ij3]we describe the results 
derived from the average [O III] A.5007 profiles, paying particular at- 
tention to how the [O IIl]A,5007 profile changes as a function of the 
parameters outlined above. In i|4] we present the results froin our 
multi-component fitting routine and calculate the fractions of local, 
optically-selected AGNs that show broad [O IIl]A.5007 lines. We 
discuss the implications of our main results in ^\ and summarise 
our findings in i|6] We adopt //y = 71 km s^' Mpc^' , £ly^ = 0.27, 
and flyv = 0.73 throughout. 



2 SAMPLE SELECTION, CLASSIFICATION AND 
ANALYSIS 

The large number of flux-calibrated galaxy and AGN spectra con- 
tained within the SDSS make it an ideal resource for deriving the 
average [O IIl]A.5007 profiles of optically selected AGNs and quan- 
tifying the proportions of the AGN sub-population that display 
broad [O IIl]A,5007 profiles. However, careful source selection is 
necessary to remove contaminants such as Hll galaxies from our 
sample of SDSS AGNs. We start this section with a description of 
how we used the results from the SDSS spectral fitting routine to 
efficiently remove those sources whose optical emission lines are 
least likely to be the result of AGN-dominated ionisation processes 
(e.g. Hll, LINERS etc.). We then describe how we used the results 
of our own emission line fitting routine to: (a) further refine our 
AGN sample, removing any remaining contaminants in the pro- 
cess and (b) measure any broad emission line components likely 
missed by the single-component fits carried out by the SDSS emis- 
sion line fitting routine. Next, we explain how we used the FIRST 
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Figure 1. Plots highlighting the line ratio and FWHMho cuts used to identify candidate Type 1 and Type 2 AGNs (shown as blue and red points, respectively; 
black points are used to indicate LINERS and Hll galaxies) in the SDSS parent sample of emission line galaxies. All values used to create both these plots 
were taken directly from the SDSS database. Left: [O IIl]^5007/Hp vs. [N Il]A6584/Ha BPT diagram of all galaxies in the parent sample. Dashed lines 
indicate the criteiia used to discriminate between candidate naiTow-line objects (i.e.. Type 2 AGNs, LINERS and HII galaxies; disciimination lines taken 
from Ho et al. 1997; Kauftmann et al. 2003). Right:[0 IIl]A5007/H(3 vs. FWHMna of the same sources. The width of the Ha line was used to discriminate 
between narrow line galaxies (i.e., candidate Type 2 AGNs, LINERS and Hll galaxies) and candidate "broad line" (Type 1) AGNs. A conservative FWHMHa 
> 600 km s^' cut was used to identify candidate Type 1 AGNs. We note that the bimodal distribution of FWHMna (a result of an apparent dearth of AGNs 
around Ha= 2000 km s"') is a feature produced by the SDSS spectral fitting routine. This bimodality disappears when the Ha profiles aie lit with multiple 
components (see i|2.3t . Contours are used in both plots to show regions of very high densities of narrow line galaxies (i.e.. Type 2 AGNs, LINERS and Hll 
galaxies). 



and NVSS radio coverage of the SDSS field to identify the radio 
loud/luminous AGNs in our SDSS sample, allowing us to explore 
the potential role of radio AGNs in driving outflows. Finally, we 
describe our spectral stacking procedure which we use as our prin- 
cipal form of analyses 

2.1 SDSS Query 

As the main focus of this study is to determine how the profile of 
the [O IIl]A,5007 emission line relates to other key AGN properties, 
we selected all extragalactic spectrosc opic sources in the SDSS- 
DR7 catalogue jAbazaiian et al.11200 9') with [O IIl]X5007 detected 
in emission (at > 3ct) by the SDSS spectral fitting routine. Since 
we rely on emission line diagnostics to identify AGNs (see i]2.2t 
we required that the Ha, HP and [N n]A.6584 lines are also de- 
tected in emission (again, at > 3o confidence). In ensuring that 
the [N Il]X6584 line is detected, sources at redshifts z > 0.4 (at 
which point [N Il]A.6584 is shifted out of the wavelength coverage 
of the SDSS spectra) are, by default, excluded from the sample. 
This query returned 239,083 SDSS targets. 

2.2 Selection of AGNs 

The vast majority of extragalactic emission line sources in the 
SDSS database are dominated by star formation rather than AGN 



activity and need to be removed from the sample prior to fur- 
ther analysis. We us ed a combination of emission line flux ratios 
("BPT" diagnostics: [Baldwin et alj[l98lh and line widths (as mea- 
sured by the SDSS spectral fitting routine) to remove the ~90 
per cent of the initial sample likely dominated by star forma- 
tion. We classified all narrow line objects (defined here as hav- 
ing FWHMHa< 600 km s"') as being probable Type 2 AGNs if 
their [O IIl]X5007/Hp and [Nll]A.6584^a flux ratios satisfy the di- 
agnostics outlined in Kauffman n etal] j2003l : see our Fig. \Q and 
identify all sources with FWHMHa> 600 km s^'as potential Type 
1 AGNs0 Our choice of FWHMea cut is a compromise between 
obtaining as complete a sample of Type 1 AGNs as possible and the 
exclusion of HlI galaxies and LINERS, whose numbers increase 
significantly at FWHMHa< 600 km s^'. Tests of a 10 per-cent 
sampling of these HlI and LINERs using our multi-component fit- 
ting routine (see i]2.3l and Appendix A) suggest a ~ 0.3 per cent 
Type-1 AGN contamination amongst these non-AGN. As such, we 
estimate we are missing ~700 Type 1 AGNs by excluding HlI and 
LINERs from further analysis. These "unidentified" Type 1 AGNs 
will be the focus of a later study. 

The above cuts identified 25,670 potential Type 1 and Type 

^ We rely on the FWHM of the Ha line, rather than Hp, to ensure we in- 
clude intermediate-class AGNs (e.g. Type 1.8, 1.9 etc.) in our Type 1 clas- 
sification 
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2 AGNs that were passed to our multi-component fitting routine 
for further analysis and possible re-classification. The remaining 
213,413 galaxies are excluded from any further analysis. 



2.3 Multi-Component Fitting Routine 

In the previous subsection we described how we used the results 
from the SDSS emission line fitting routine to crudely select sam- 
ples of Type 1 and Type 2 AGNs in the SDSS database. However, 
the single gaussian (i.e., one per emission line) fits performed by 
the SDSS pipeline cannot account for the complex emission line 
profiles that are common among AGN spectra. Instead, we devel- 
oped our own multi-component fitting routine to accurately mea- 
sure these line profiles. We used this routine to re-measure the pro- 
files of the Ha, HP, [N Il]A,6548,6584 and [O IIl]A,5007 emission 
lines in all 25,670 spectra identified in ^2.2\ The results from this 
multi-component fitting routine were used to (a) obtain more robust 
AGN classifications and (b) identify any broadening or shifting of 
the [O IIl]A.5007 emission line. A detailed description of our multi- 
component fitting routine is given in appendix A. To summarise, 
up to two gaussians are fit to each of the forbidden [O IIl]A,4959, 
[O IIl]X5007, [N Il]X6548 and [N Il]A,6584 emission fines (to ac- 
count for any line asymmetries and/or broad components) and up 
to three gaussians to each of the permitted Ha and Hp lines (up to 
two gaussians to account for the narrow component of these per- 
mitted lines and one gaussian to account for the broad component 
in the case of type I AGNs). The two gaussians used to model the 
narrow component of these permitted lines are fixed to have the 
same profile as the forbidden [O III] and [N II] lines. 



2.4 Source Classification 

We used the results from our multi-component line fitting routine 
to refine the classifications of the AGNs in our sample beyond that 
which is capable using the single component fits performed by the 
standard SDSS pipeline. In what follows, we outline how we sep- 
arated our sample into broad and narrow line Seyfert Is (hereafter, 
BLSl, NLSl respectively, and referred to collectively as Type 1 
AGNs), Type 2 AGNs and those sources that, following analysis 
with our multi-component fitting routine, were identified as being 
non-AGN dominated (and excluded from further consideration). 



2.4. 1 Type 1 AGNs (BLSls and NLSls) 

For each source we used the FWHM of the broadest Ha com- 
ponent to identify Type 1 AGNs and distinguish between BLSls 
and NLSls. By using Ha rather than HP as the discriminator, 
we will avoid cl assifying intermediate type AGNs (e.g., 1.8, 1.9; 
IOsterbrocklll98 l') as NLSls since, by definition, intermediate types 
show evidence of a broad component in the Ha line (but not neces- 
sarily in the HP line). A source was classified as a Type 1 AGN if an 
extra gaussian (i.e., beyond that needed to fit the forbidden lines) 
provided a significantly better fit (at >99 per cent confidence) to 
the Ha line, and that gaussian constitutes at least 50 per cent of the 
total emission fine flux and has FWHMHa> 600 km s"' (10,548 
sources). The > 50 per cent cut excludes 995 AGNs that our fit- 
ting routine measures FWHMho(> 600 km s^', 717 (i.e., ~72 per 
cent) of which are subsequently classified as type 2 AGNs. These 
995 AGNs represent < 10 per cent of either Type Is or Type 2s 



and, subsequently, their reclassification has little impact on our re- 
sults yet ensures that we only include genuine type 1 AGNs in that 
group. 

As per the classical, yet arbitrary, definition of iGoodriclil 

( Il989h . all of those Type 1 AGNs whose broadest Ha compo- 
nent has a FWHMHa> 2000 km s^' were sub-classified as BLSls 
(9,455 sources). The remaining Type 1 AGNs (1,093 sources) were 
classed as NLSls (of which 1,024, or ~93 per cent, have broad Ha 
components that satisfy 1000 km s"' <FWHMHa<2000 km s"'). 

2.4.2 Type 2 AGNs 

We used the lKauffmann et alj j2003h emission line (BPT) diagnos- 
tics to identify Type 2 AGNs in the remaining sample (i.e., those 
that have not already been identified as Type 1 AGNs). Where two 
gaussians are fit to the [O IIl]X4959,5007 and [N n]A,6548,6584 
lines, we use the combined flux of these components. In such cases, 
we use the combined flux of the two narrow components of the HP 
and Ha lines. Similarly, when only one gaussian is needed to fit 
the [O IIl]X4959,5007 and [N Il]X6548,6584 fines, we only use the 
flux of the single narrowest gaussian of the Ha and HP lines (i.e., 
if present, we always exclude any Ha and Hp broad components 
when calculating their narrow-line flux). In this way, we classify 
13,716 Type 2 AGNs. 

2.4.3 Non-AGN Dominated 

Any narrow line objects (i.e., non-Type 1 AGNs) which, after being 
analysed by our multi-component fitting routine, were not classified 
as Type 2 AGNs were removed from the sample (1,406 sources), 
leaving 24,264 confirmed AGNs (of either Type 1 or Type 2). Af- 
ter analyses with our emission-line fitting routine, these rejected 
sources were found to predominantly lie in the LINER region of 
the BPT diagram and have FWHMHa> 600 km s^' as measured 
by the SDSS pipeline. 

2.5 1.4 GHz Radio Coverage 

One of the key aims of this study is to determine how the kinematics 
of the [O III] A.5007 emitting gas relates to the radio properties of the 
AGN. To quantify the radio properties of our sample we performed 
a three-way cross matching between our optical (SDSS) catalogue 
and the NVSS and F IRST radio catalogues, largely following the 
procedure outlined in lBest et al. I (2005.) . As such, we only use radio 
fluxes from the NVSS catalogue, but use the higher-spatial resolu- 
tion information provided by the FIRST survey to ensure a low 
level of con tamination within our optical-radio matched catalogue. 
Contrary to iBest et al. I( l2005h . the goal of this study is not to ob- 
tain a highly complete sample of radio sources, but instead to mea- 
sure the radio properties of optically selected AGNs. As such, we 
considered all NVSS sources with 1.4 GHz fluxes > 2 mjy; this 
lower limit is a factor of 2.5 below that used by Best et al. (200^. 
However, we note that the integrated 1.4 GHz flux densities of all 
the sources in our catalogue have signal-to-noise ratios > 3. For 
those sources not detected by the NVSS, we assumed an integrated 
1 .4 GHz flux density upper limit of 2 mJy. 

Of the 10,548 Type 1 and 13,716 Type 2 AGNs in our SDSS 
sample, 535 (i.e., ~5 per cent) and 1453 (i.e., ~1I per cent), re- 
spectively, are found to have radio counterparts (representing ~8 
per cent of all AGNs in our sample). The discrepancy between the 
fraction of Type I AGNs with radio detections compared to that 
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of Type 2 is likely to be due to the redshift distributions of these 
classes of galaxies in our sample. Since they are brighter at optical 
wavelengths, Type 1 AGNs are more likely to satisfy the appar- 
ent magnitude cut for spectral follow-up by the SDSS. The higher 
average redshift of the Type 1 AGNs means a smaller fraction of 
them are detected in the NVSS radio survey (which has roughly a 
constant flux limit). 

We assume a radio spectral index of a = 0.8 when calcu- 
lating the rest-fr ame 1.4 GHz radio luminosity, where Fy v^" 
dlbar et al.ll201(]t) . 

2.6 Spectral stacks and average [O iii]A.5007 profiles 

In addition to our fits to individual spectra we also rely heavily 
on spectral stacking to derive the average [O IIl]X5007 profiles of 
optically-selected AGNs. By binning the AGN sample in terms of 
luminosity, Eddington ratio and radio properties we can use these 
stacks to isolate general trends in the [O IIl]A.5007 profile with re- 
spect to these parameters. The benefit of spectral stacking is that it 
is less susceptible to systematic biases than results derived solely 
from individual fits (e.g., the low S/N of fainter AGNs affecting 
our ability to detect broadened line components) and, by consider- 
ing relative changes in the average profile of the [O IIl]A.5007 line, 
to measures of the systematic velocity of the AGN. In S|4]we use the 
results from our spectral stacks to motivate further analyses using 
the individual AGN spectra. 

The spectral stacks were built by coadding the continuum- 
subtracted spectra provided by the SDSS database. Each spectrum 
was shifted into the rest frame using the [O IIl]A.5007 redshift de- 
rived from our fits (assuming that the narrowest [O IIl]A.5007 com- 
ponent lies at the systematic redshift). We note that the redshift de- 
rived from the [O IIl]A,5007 line is in very good agreement with 
the redshift provided by the SDSS database - a gaussian fit to 
the Az = ZSDSS ^ Z[o III] distribution reveals an average offset of 
Az= (0.9 ±2.1) X 10"'* (i.e., 26 ±63 km s"'); indeed, our tests 
show there is no significant impact in our results if we adopt the 
SDSS-reported redshifts instead. We take a bootstrapping approach 
to calculate the error on the average emission line profiles. We se- 
lect, at random, a third of the sample in each of our bins and stack 
their spectra. We repeat this process 1000 times and calculate the 
standard deviation at each wavelength of the resulting distribution 
of stacked spectra. This standard deviation is then adopted as the 
error on the stacked line profiles. 

When stacking as a function of radio properties (i.e., Lj 4 qhz 
and -RRad)i we only include radio-undetected sources when (a) their 
upper limits are consistent with the upper boundary of the bin 
and (b) there is no lower boundary on the bin (e.g., -Li,4ghz< 
10^2 W Hz- 1 andiJRad< 10"'). 



3 RESULTS FROM SPECTRAL STACKING 

In our analyses we explore the kinematics of the AGNs' narrow line 
regions as portrayed by the profiles of the forbidden [O IIl]X5007 
emission lines. Using stacked spectra, we search for general trends 
in the profiles of the [O IIl]X5007 lines with respect to various key 
AGN parameters (specifically, bolometric luminosity, Eddington 
ratio, radio loudness, radio luminosity and AGN orientation). In ij4] 
we use the results from the emission line fitting routine described 
in 323] (see also appendix [At to estimate the proportions of opti- 
cally selected AGNs that have either very broad or highly shifted 
[O IIl]X5007 lines or components. 
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Figure 2. Average [O ni]X5007 profile of Type 1 AGNs in our sam- 
ple derived from spectral stacking. This profile shows a strong blue wing 
(bluewai'd asymmetiy) that is well modelled by a broad (i.e., FWHM= 
851 km s"') second gaussian component that is offset from the nan'ower 
component (i.e., FWHM= 335 fan s"') by 148 km s"'. The fit produced 
by combining these two components is shown by the dashed line. 

We initially focus our analyses on Type 1 AGNs. Assuming 
outflows are driven radially outwards perpendicular to the plane of 
the d usty torus/accreti on disk, then from AGN unification models 
(e.g., lAntonuccll 1993b any outflowing material driven by the AGN 
(either by jets or radiatively driven winds) should show a stronger 
blueshift (asymmetry) in the spectra of Type I AGNs (compared 
to Type 2s) as the outflow will be directed more along our line 
of sight. As such, if the AGNs are indeed driving outflows, and 
these outflows emit [O IIl]X5007, we should expect Type 1 AGNs 
to show stronger [O IIl]X5007 line shifts. In i]3.3l we compare the 
results derived from our Type 1 population with the Type 2 AGNs 
in our sample to explore the influence of AGN orientation on the 
observed emission line profiles and measured gas kinematics. 

3.1 The average [O Iii]l5007 profile of Type 1 AGNs 

In fig. [2] we present the average [O IIl]A.5007 profile derived from 
stacking all Type 1 AGNs in our sample. This profile displays 
a prominent blue wing and is well fit using two gaussian com- 
ponents; a narrow (i.e., FWHM= 335 km s"') component and a 
broader (i.e., FWHM= 851 km s"') component that is blueshifted 
by 148 km s"' with respect to the narrow component. The flux con- 
tained within the broader component represents ~45 per cent of the 
total line flux. 

In the following subsection we consider whether the profile 
of the [O III] ^5007 line changes as a function of bolometric lumi- 
nosity (Lagn)' Eddington ratio (A.Edd)> radio loudness (-/?Rad)> radio 
luminosity (Li.4ghz) and orientation. We do this by binning the 
sample in terms of these properties and use spectral stacking to de- 
rive the average [O IIl]X5007 emission line profiles of the AGNs in 
each bin. 
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Figure 3. Average [O IIl]^5007 profiles of Type 1 AGNs binned according to (a) [O III] luminosity, (b) Eddington ratio, (c) radio loudness and (d) radio 
luminosity. The grey curve in plot (b) shows the average profile of the highest Eddington ratio AGNs after excluding those spectra that show prominent Fe II 
emission that can contaminate [O IIl]X5007. For presentation purposes we do not show the errors on the stacked profiles which, with the exception of the 
highest Eddington ratio bin in panel (b) (see representative error bar), are smaller that the points. 



3.2 Correlations between the average [O iii]X5007 profile and 
key AGN properties 

3.2.1 [Om\k.50Ql luminosity 

If outflows from AGNs are driven by radiative forces, and these out- 
flows manifest themselves as blueshifted [O IIl]A,5007 wings, then 
we may expect the observed properties of these wings to change 
as a function of the AGN luminosity. The [O IIl]A.5007 luminosity 
(hereafter, m]) is thought t o be a proxy measure of the bolo- 
metric AGN lum inosity (e.g. Maiolino & Riek jl995l : lRisaliti et al.l 
1 19991 J Bassani et al. 1999 ; Heckman et al. 2005), so here we inves- 
tigate the relationship between [O IIl]X5007 profile and lu] • 

In the top-left panel of fig. [3] we present the average 



[O IIl]X5007 profiles of Type 1 AGNs in our sample divided into 
bins of log(Ljoin])5 with bin intervals of Alog(Ljo in])=l over 
LjO iu]= 10^" ^ 10"*^ srgs s^'. We calculate Ljq hi] using the re- 
sults from our multi-component fitting routine, assuming the red- 
shift calculated from the [O IIl]A,5007 line; when two [O ni]A.5007 
components are present, Ljq m] is the combined luminosity of both 
components, although we note that our results remain unchanged 
when only the luminosity of the narrow [O IIl]A,5007 component is 
used. We use the Balmer decrement (using the ratio of the narrow 
Hg and Hp compon ents calculated from our fits) and assume the 
ICardelli et al] ( Il989l) reddening curve to correct for the effects of 
reddening on Ljq m] ■ 

The stacked [O IIl]A,5007 emission line profiles of Type 1 
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Figure 4. Average [O IIl]^5007 profiles of Type 1 AGNs, accounting for the correlation between i.|o iii| and Li.4gHz in an attempt to isolate which has 
the greater affect on the line profiles. Panels (a) and (b) show the results of stacking as a function of Lpm], while keeping Li.4gHz roughly fixed (to 
< 10-^ W Hz"' and Kfi'^ W Hz"', respectively). Alternatively, panels (c) and (d) show the results of stacking as a function of Li 4 ghz keeping Lp m] roughly 
fixed (to ~ 5 X lO**' ergs s"' and ~ 2 x lO**" ergs s"', respectively). The greatest change in the profiles of the [O IIl]X5007 line is seen when Li 4 ghz varied. 
As with fig.[3]we show representative error bars when they are larger than the plotted points (see panels (c) and (d). 



AGNs display increasingly prominent blue wings with increasing 



-[O HI] 



However, before interpreting this result as evidence that 



radiative forces are largely responsible for the kinematics of the 
[O IIl]A.5007 emitting gas, it is important to account for the known 
positive correlation betwe en Lp jitj and A GN ra dio power (e.g., 
Baum & Heckman' 'l989|^ iRawlings et li\ 1 19891: IZi rbel & Baunj 
1995l : lTad hunter et al. 199ll: lwills et al.l2004l:lde Vries et al. 2007). 
For this, we perform the same experiment as described above (i.e., 
stacking as a function of Ljg m] ) however this time controlling for 
radio luminosity to ensure that the median radio luminosities of 
the stacks are consistent to within a factor of 3 of each other (see 
Table [T). We have also attempted to keep the number of sources 



in each stack roughly the same by randomly selecting, without re- 
placement, sources from each bin. The numbers in each stack are 
therefore constrained by two effects: (i) the number of sources in 
the bin containing the smallest number of AGNs and (ii) the ability 
to select sources from each bin with radio luminosities that are con- 
sistent to within a factor of 3. The combination of these two effects 
is why the numbers of stacked sources are not exactly the same in 
each of the bins. 

The results of this experiment are shown in the top two pan- 
els of fig. |4] Here, we have divided the sample into radio weak 
(i.e., Li 4 GHz< 10^^ W Hz"', including any upper limits that also 
satisfy this criterion: panel a) and radio powerful (i.e., Li 4ghz> 
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10^3 W Hz-i, excluding all upper limits; panel b). When the 
stacks are controlled for radio luminosity, the change in average 
[O IIl]X5007 profile as a function of L[o ni] is much less pronounced 
in both the radio weak and radio strong samples|3 We do note, how- 
ever, that a slight increase in the prominence of the [O IIl]A.5007 
blue wings with increasing Ljq m] still remains in both the radio 
strong and weak subsamples. However, comparing the top two pan- 
els of fig. |4] reveals a striking difference in the average profiles 
of radio weak and radio strong AGNs, with all the radio strong 
stacks having broader [O IIl]A,5007 profiles than their equivalent 
radio weak stack. We explore this further in i]3.2.3l 



3.2.2 Eddington ratio 

Previous studies of the [O IIl]X5007 line profile have concluded 
that the Eddington ratio (i.e., "k^a = iAGN/^Edd> where L^^^ 
is the Eddington luminosity of the central SMBH) of the AGN 
plays a key role in dictating both the FWHM and blueshift of the 
[O II ljA.5007 line, inc luding whether it displays a strong blue wing 
fe.g.. lGreene & Hol2005a) . Here, we consider whether our analyses 
confirm that the [O IIl]X5007 profile is correlated with the Edding- 
ton ratio of the AGN. 

To estimate the Eddington luminosity of the AGN we use the 
SMBH mass calculated fro m the FWHM and lu minosity of the Ha 
emission line (see eqn. 6 of lGreene & Holl2005b) . We only calcu- 
late LEdd for Type 1 AGNs and use only the FWHM of the broadest 
Ha gaussian. Here, the bolometric luminosity (i.e., Lagn) is es- 
timated from the luminosity of the Ha line using a combination 
of the L[ia-Ls\nn and L<,\o f)-LAc,N correction factors des cribed in 
iGreene & Horj2005ah and iNetzer & Trakhtenbro3 (l2007h . respec- 
tively. These combine to give an expression for the Eddington ratio 
of: 



Table 1. Sample properties after accounting for the correlation between 
^|o m] and Li,4 ghz- 



'>-Edd = 



i-AGN 



^Edd 



'0.43/10 



10"*^ ergs s 



0.31 



FWHMHa 

10^ kms-1 



-2 



(1) 



where f\Q is the dimensionless L^\o(\-L agn bolometric coiTection 
factor (see lNetzer & Trakhtenbrotll2007b in units of 10. We assume 
/lO = 1 throughout, although note that this has no effect on our 
final results as we are only concerned with relative changes in the 
[O IIl]X5007 profile. 

The average profiles of the Type 1 AGNs in our sample, 
stacked in terms of "k^id ars presented in the top-right panel of fig. 
[3] We find no significant difference in the average [O ni]A.5007 
profiles over the range XEdd= 3 x 10^"' — 3 x 10^^, while the 
[O IIl]A,5007 profiles of A,Edd= 3 x 10"^ _ i ^GNs display in- 
creasingly prominent blue wings with increasing k^a. The aver- 
age spectrum of AGNs in the highest "k^a bin that we consider 
(i.e., A.E(id= 0.3 — 1) shows evidence of contamination from Fe 
11 (beyond, e.g., ~ 5020A) that could alter the apparent profile 
of the [O lll]A.5007 line0 When we remove those AGNs in the 
highest Xeaa bin that show evidence of having strong Fe II emis- 
sion at other wavelengths in their spectrum, the resulting stacked 
[O lll]A.5007 profile remains largely unchanged. As with our stacks 



^ We note that we do not calculate the average profile for the highest Lp m] 
bin (i.e., Lp m]= lO^^ — 10*'' ergs s^') of the radio weak AGNs, nor the 
lowest Z,(oiii| bin (i.e., L[oiii]= lO'*" — Itf*' ergs s^') of the radio strong 
AGNs, as there are too few sources to provide reliable averages - a conse- 
quence of the correlation between Lp III] ™d radio power 

There is a known c orrelation between the eq u ivalent widths o f the Fe II 
complex and Xedd, see lBoroson & Greenl l ll992l) ; lBorosonl ( l2002h . 



i-Rad 


^[OIII] 


N 


Med. LRad 


Med. Lpuii 






\i) 






< 10 


10 — 10 


129 


< 23 


40.8 


< lO^"" 


10^' - 10*^ 


129 


< 23 


41.6 


^ in23 


1 n4^ 1 n43 
10 — 10 


129 


< 23 


42.3 


> 10^3 


1041 _ 1042 


64 


24.1 


41.7 


> 10^3 


10« _ 1043 


77 


24.0 


42.5 


> 10^3 


1043 _ 1044 


77 


24.1 


43.3 


<1022 


<10''2 


33 


< 22 


41.5 


1022 _ 1023 


<1042 


33 


22.6 


41.7 


1023 - 102^* 


<1042 


32 


23.4 


41.8 


102-* - 1025 


< 1042 


33 


24.3 


41.7 


<1022 


>1042 


39 


< 22 


42.1 


1022 _ io23 


>1042 


42 


22.7 


42.3 


1023 _ 1024 


>1042 


42 


23.5 


42.4 


1024 _ 1025 


>1042 


42 


24.3 


42.4 


1025 _ io26 


>1042 


37 


25.5 


42.4 



Notes: The first six rows are the samples that have been selected to have 
roughly constant Li .4 gHz but increasing Lp m] ■ The lower nine rows are 
those samples that have been selected such that Lp m] remains roughly con- 
stant, but Li 4 GHz increases. (1) Range of L[ 4 ghz: (2) range of Lp iii|; (3) 
the number of sources in each subsample: (4) the median radio luminosity; 
(5) the median [O III];>l5007. 

as a function of Lp mj , we must use caution before directly as- 
sociating this broadening with increased A-Edd. as there may be 
other, more fundamental, drivers. To test for this, we performed a 
similar experiment as outlined at the end of i]3.2.1l splitting the 
high A.Edd AGNs into bins of L[o m] and Li.4GHzWhile control- 
ling for A-Edd- Unfortunately, there are too few AGNs in our high- 
est A-Edd bin (i.e., A,Edd> 0.3) to perform a meaningful experiment 
(there are only six Type 1 AGNs in our sample with XEdd> 0.3 and 
i'l.4GHz> 10^"' W Hz^'), so we resort to A.Edd> 0.1). The results 
of this experiment is shown in fig.|5] As with the full Type 1 AGN 
sample, we see a broadening of the [O lll]A.5007 line with Lp mj 
for k^dd> 0. 1 AGNs, but note that these stacks have not been con- 
trolled for Li 4 GHz> so some of this broadening may be related to 
the correlation between Lp mj and 4 ghz- Indeed, we see a strik- 
ing difference between the [O lll]A.5007 profiles of high and low 
Li 4 GHz AGNs with high k^^d^ suggesting that even in these high 
A-Edd AGNs, mechanical processes are playing an important role in 
the disturbing the [O lIl]X5007 emitting gas. However, we stress 
that this only applies to A,Edd> 0.1 AGNs, as the small numbers of 
A,Edd> 0.3 AGNs in our sample means we cannot perform this ex- 
periment in the highest k^dd bin where we see the strongest blue 
wing. 



3.2.3 Radio loudness and luminosity 

In addition to radiative forces, it is also thought that radio jets (ei- 
ther compact or extended) may play a significant role in disturbing 
the interstellar medium arou nd an AG N, potentially inducin g feed- 
back (e.g.. Whittle 1992; Nelson & W hittle 1996; Tadhunt er et al.l 
l200ll ; lHolt et ah.2003i , .200^ . Indeed. iNesvadba et al . |200§') have 
suggested that the spatially resolved broadening of the [O 111] 15007 
lines in the spectra of high redshift radio galaxies are the result of 
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Figure 5. Average [O IIl]^5007 profiles of high Eddington-ratio (^Edd> 0.1) Type 1 AGNs, split as a function of [O IIl]X5007 luminosity (left) and radio 
luminosity (right). Similar to Fig.|4] we have stacked in different Lp ni] ™d Li 4 ghz bins (see keys), matching ^Edd in sach bin to keep the average 
roughly the same between the two bins. As with fig.|3]we show representative en'or bars when they are larger than the plotted points. 



interactions between a radio jet and the NLR. If this the case, then 
we may expect the profiles of the [O ni]A.5007 line to be related to 
the radio properties of the AGNs. 

In the lower two panels of fig. [3] we present the average 
[O IIl]X5007 profiles of the Type 1 AGNs in our sample, stacked 
according to radio loudness (defined here as i^Rad = ^1.4 GHz/^b; 
where Fi 4 ghz and are the flux densities measured in the 
1.4 GHz band and the rest-frame optical B-band, respectively) 
and radio luminosity (at 1.4 GHz; hereafter, Li.4ghz)Ij The av- 
erage [O IIl]A.5007 profiles show evidence of only a modest change 
across over three orders of magnitude in i^Rad (i-e., ^Rad< 10 to 
10^ - 10"^). As an aside, we find that the average [O IIl]A.5007 
lines of radio detected (iiTespective of loudness) AGNs have pro- 
files similar to those AGNs with the highest L|o m] values (i.e., 
LjO in| > 10^^ srgs s^ ' ). However, this result is not surprising when 
we consider that the radio detected AGNs in our sample typi- 
cally have high [O IIl]A.5007 luminosities. For example, 75 per 
cent (i.e., 654/871) of radio detected Type 1 AGNs have Ljoiii|> 
10^^ ergs s^'). 

We also investigated whether any trend exists between the av- 
erage [O IIl]A.5007 profile and radio luminosity (as opposed to ra- 
dio loudness). From our spectral stacks we find that there is a non- 
monotonic change in the average [O IIl]A.5007 profiles over the four 
orders-of-magnitude range in Li 4 ghz spanned by our sample; the 
it. 4 GHz= 10^^ - 10^4 W Hz- 1 andLi,4 ghz= lO^'* - 10^5 WHz-' 
stacks display a marginally stronger blue wing compared to the 
it.4 GHz= 10^^ - 10^3 WHz-i andLi.4 ghz= 10^^ - 10^^ Wnz"' 
stacks. 

As outlined in i)3.2.1l it is important that we consider the in- 
fluence that the correlation between Lp m] and L1.4 ghz rnay have 

^ B-band fluxes were calculated by passing the rest-frame SDSS spectra 
through a synthetic Buser B-band response filter 



on our results. For this, we employ a similar method as that used 
in i]3.2.1| This time, however, we control for L|o m] while stack- 
ing as a function of L1.4 ghz> ensuring that the median L[o in| of 
the stacked sources are consistent to within a factor of 3. The 
results of this experiment are shown in the lower two panels of 
Fig. |4] Here, we have separated the sample into high and low 
[O IIl]A.5007 luminosity subsamples (i.e., Z'[oiii]> 10'*^ ergs s^' 
and L[oiii]< 10^^ ergs s^', respectively). In both subsamples we 
see a clear change in the average [O IIl]A,5007 line profiles as a 
function of Li.4ghz- We note again, however, that this increase 
is not monotonic. Instead, the [O IIl]A.5007 profile broadens with 
increasing 4 ghz from Li.4ghz< 10"^ W Hz^' to Li.4ghz= 
10^3 - 10^"^ W Hz^i (L1.4 GHz= 10^^ - 10^5 W Hz^i in the case 
of L[o ni]> IC*^ ergs s^' AGNs), but then decreases in width as 
Li 4 GHz increases further (to Li 4 ghz= 10^"* ^ 10^^ W Hz^' and 
il.4 GHz= 10^^ - 10^'' W Hz^i for the m] < 10"^^ ergs s"' and 
Ljo iu]> 10^^ ergs s^' subsamples, respectively). 

3.2.4 The average [O IIl]A,5007 profiles of Narrow -line Seyfert Is 

It has previously been suggested that Type 1 AGNs with nar- 
row Balmer lines (i.e., < 4000 km s^') and, in particular, NLSls 
show a tendency to have more prominen t [G IIl] A,5007 blue 
wings or blueshifts c ompared to BLSls (e.g., lMarzian i et al. 200^; 
iKornossa et alj|2008l ) - a feature previously attributed to the higher 
Eddington ratios of NLSls. We can test for this using our sample of 
Type 1 AGNs. In fig.|6]we present the stacked [O IIl]A.5007 profile 
of the NLS Is in our sample compared against that of BLS Is. Here, 
we have stacked the same number of randomly selected BLSls 
as there are NLSls in our sample (i.e., 1096). We find no signifi- 
cant difference between the average [O IIl]A.5007 profiles of BLSls 
and NLSls in our sample. As a check, we also show in this plot 
the stacked Ha profiles, confirming that our sample of NLSls do. 
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Figure 6. Average [O IIl]^5007 profile (left) and Ha profile (right) of the BLSl and NLSl AGNs in our sample. Inset we show the distribution of Eddington 
ratios for the NLSls (black histogram) and the BLSls (blue histogram) in our sample, indicating that the NLSls have systematically higher Eddington ratios. 
We find no evidence to suggest that, on average, NLSls have [O in]X5007 lines that are either broader or more blueshifted than their BLSl counterparts. 



indeed, have significantly narrower permitted emission lines (i.e., 
broad component FWHMea of 4460 km s^'and 1590 km s^' ; typ- 
ical for these respective classes of object). Furthermore, the NLSls 
in our sample have systematically higher Eddington ratios than 
their BLS 1 counterparts, as highlighted in the inset plot showing 
the distribution of A-Edd for the BLSl and NLSl samples. Finally, 
we note that the profile of the stacked [O ni]A.5007 is largely inde- 
pendent of our criteria for identifying NLSls; there is no significant 
change in the [O IIl]A.5007 profile when we change our NLSl se- 
lection criteria to FWHMHa< 1500 km s"' or < 2500 km s"'; 
in other words the average [O ni]A.5007 profile is not dominated 
by NLS 1 AGNs with FWHMea close to the (arbitrary) boundary 
between NLSls and BLSls. 



3.3 Comparison of [O ill] A,5007 profiles of Type 1 & Type 2 
AGNs 

Up to this point, we have only considered the [O ni]A.5007 pro- 
files of Type 1 AGNs. The observed differences between Type 1 
and Type 2 AGNs are now generally regarded as being the result 
of obscuration caused by gas and dust alo ng our line of sight to 
their central regions (e.g., lAntonucdlll993l) . It is therefore of inter- 
est to consider Type 2 AGNs in order to gain insight into the role 
that AGN orientation and obscuration plays in the profiles of the 
[O IIl]A,5007 emission lines. 

In fig. |7]we present the average [O ni]A.5007 profile derived 
from stacking all 13,713 Type 2 AGNs in our sample. For com- 
parison, we also include the average [O ni]A.5007 profile of the 
Type 1 AGNs in our sample (presented previously in fig. H). Com- 
pared to the average [O IIl]A,5007 profile of the Type 1 AGNs, the 
Type 2s display a less prominent blue wing, (i.e., is more sym- 
metrical). Since the [O IIl]A,5007 line broadens with increasing 



LjO III] (but see i]3.2.3t . the narrower Type 2 profile may in part 
be the result of the lower average Ljq m] of the Type 2 AGNs in 
our sample compared to Type Is (due to dilution from the strong 
continuum in the latter; median -L[oiii]= 6 x 10"*' ergs s^'and 
= 1.8 X 10^^ ergs s^', respectively). Having said that, as with the 
stacks of the Type 1 AGNs, the [O III] A,5007 profiles Type 2 AGNs 
become increasingly broad with increasing Ljq m] ^ but remain more 
sy mmetric than the Type 1 stacked profiles, as recently reported 
bv lVaona et all fa 12) (see fig.l] left-hand panel). Similarly, we 
note that the [O IIl]A.5007 profiles of the Type 2 AGNs follow the 
roughly the same trend with L\n ghz as the Type Is (i.e., increas- 
ing in width to Li.4ghz= 10^"' — 10^^ W Hz^', then narrowing 
to higher Lj 4 ghz) while, again, remaining more symmetrical (see 
fig. [8] right-hand panel). The more symmetrical [O IIl]A.5007 pro- 
files of Type 2 AGNs are expected if the central engine is driv- 
ing an [O IIl]A,5007 emitting outflow along the axis perpendicular 
to the plane of the obscuring torus. In this model, orientation ef- 
fects mean that outflows from type 1 (i.e., unobscured) AGNs will 
appear more strongly blueshifted than those from type 2 (i.e ., ob- 
scured) AGNs. This is consistent with lHarrison et ail J2012h . who 
reached the same conclusion using spatially resolved (i.e., integral 
field unit) spectroscopy. 



4 RESULTS FROM INDIVIDUAL EMISSION LINE FITS 

We have shown that the average [O IIl]X5007 emission line profile 
of optically selected Type 1 AGNs displays a strong blue wing and 
that, once the correlation between L[o m] and Li.4ghz has been 
taken into account, the width of the average [O IIl]A.5007 is most 
strongly related (though not in a monotonic sense) to the radio lu- 
minosity of the AGN. We see no evidence that the average profile 
of the [O IIl]A.5007 line changes strongly as a function of radio 
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Figure 8. Average [O IIl]^5007 profiles of Type 2 AGNs binned according to [O in]^5007 luminosity. As is the case for Type 1 AGNs (see fig.|2}, we find 
that the average [O IIl]^5007 profile of Type 2 AGNs increases in width with increasing [O in]X5007 luminosity. However, the Type 2 profiles show much less 
prominent blueward asymmetries compared to the average [O IIl]X5007 profiles of Type 1 AGNs. As with fig.|3]we show representative error bars when they 
are larger than the plotted points. 
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Figure 7. Average [O III] ^5007 profile of the Type 2 AGNs in our sample, 
derived from spectral stacking. This average profile shows less prominent 
blue wings compared to the average Type 1 profile (shown in grey). The 
absorption either side of the [O IIl]^5007 line is attributable to the stellar 
features in the host galaxy component which are stronger in Type 2 AGNs 
due to the central engine continuum being obscured from view. 



loudness or Eddington ratio, and find that ni] (assumed to be a 
proxy for Lagn) has a relatively minor impact on the average pro- 
file (compared to the radio luminosity). There is no difference in 
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Figure 9. The fraction of AGNs in our samples with FWHMavs greater than 
a given value (shown on the abscissa). [O IIl]X50()7 lines with FWHMAvg 
> 1000 km s^' are Ri5 times more prevalent among AGNs with Li 4 gHz> 
10^' W Hz^' compared to the overall AGN population. We have included 
upper limits in the Li 4 gHz< 10^' W Hz"' sample when they are consistent 
with this limit. 



the average [O IIl]X5007 profiles of BLSls and NLSls, while Type 
2 AGNs have, on average, more symmetric [O IIl]A,5007 emission 
lines compared to Type Is but their widths follow the same general 
trends. 

In this section, we further explore the results derived from the 
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spectral stacks by using the results from our fits to the [O IIl]A.5007 
emission line of the individual spectra in our AGN sample. In doing 
so, we will demonstrate that the vast majority of all AGNs have rel- 
atively narrow average [O IIl]A.5007 profiles, while it is the AGNs 
with lO^-' <Li.4GHz/W Hz"^< 10^^ that dominate the numbers 
with the broadest [O IIl]X5007 lines. 

We use the results from our multi-component fitting routine to 
quantify the proportions of optically selected AGNs with broad- 
ened [O IIl]A.5007 lines. Our fitting routine fits up to two gaus- 
sian components to the [O IIl]X5007 lines. However, rather than 
focus on the properties of single components, we calculate the flux- 
weighted average FWHM of the [O IIl]>.5007 line: 

FWHM^,,^ = ((FWHMA/k)^ + (FWHMbFb)^)! (2) 

where Fj^ and are the fractional fluxes contained within the two 
fitted components, A and B. By primarily focussing on average 
FWHMs, we avoid arbitrary definitions such as the threshold above 
which a broad component constitutes an important contribution to 
the overall flux of the emission line, or the point beyond which 
a component is considered "broad". This characterisation of the 
width of the [O IIl]X5007 line also ensures that all AGNs, whether 
their [O III] ^5007 lines fitted either with one or two gaussians, 
can be compared on an equal footing (i.e., including, for example, 
AGNs in which the [O IIl]A.5007 line is broad, but satisfactorily fit 
with a single gaussian). 

By measuring the [O IIl]X5007 profiles of individual AGNs we 
can quantify our sample in terms FWRM^vg- In fig.|9]we show the 
fractions of our sample with FWHMAvg above a given value. It is 
clear that AGNs with relatively narrow [O III] A,5007 lines dominate 
our sample, with the two-thirds of all AGNs having FWHMAvg < 
400 km s"'. The median FWHMAvg of L[oiu]> 10"*^ ergs s"' 
AGNs is marginally higher than that of our whole AGN sample 
(~ 340 km s^' and ~ 330 km s^'), although this may in part be 
due to a broad component being easier to identify among the higher 
S/N spectra of these high L|o ni| AGNs. Approximately 17 per cent 
of all AGNs in our sample have FWHMAvg > 500 km s^' , which 
is roughly the lower FWHM limit adopted by previous studies of 
AGN outflows when selecting targets for follow-up resolved (i.e., 
long-slit or i ntegral-field) spectro scopic studies of AGN outflows 
(see Fig. 1 of lHarrison et al.ll2012l) . Also shown in fig.|9]is the cu- 
mulative fraction for Li 4 ghz> 10^"^ W Hz^' AGNs, which indi- 
cates that the incidence of extremely broad [O IIl]A,5007 lines (i.e., 
FWHMav»> 1000 km s"') is 7»5 times higher among these ra- 
dio luminous AGNs than among AGNs in general (although the 
fractions are still small; e.g., ~ 5.5% and ~ 1.1, respectively, have 
FWHMAvg > 1000 km s-'). 

In fig. [To] we plot the distribution of FWHMAvg as a function 
of L|o in] . In the leftmost panel we plot all AGNs in our sample 
(including both type 1 and type 2 AGNs). As expected from the 
cumulative fractions (fig.|9}, AGNs with FWHMAvg< 500 km s"' 
dominate this plot at all L[o values. The upper envelope of this 
distribution increases with L|o m] up to L|o in]~ 10^^ ergs s^', al- 
though we reiterate that this may partly be explained by the low 
S/N ratios of low luminosity AGNs' spectra making it difficult to 
detect a broad emission line component. Another possibility high- 
lighted by the results of our spectral stacks is that this increase is 
driven by the L|o ni|~^l.4 GHz correlation coupled with an increase 
in FWHMAvg with L1.4 ghz- To test for this, we have split the sam- 
ple according to Li 4 ghz in the other two panels of fig. llOKshow- 
ing < 10^3 W Hz^l and > lO^^ W Hz^^n the middle and right- 
hand panels, respectively). We include AGNs with upper limits on 



^1.4 GHz (provided they are < 10 W Hz^ ) in the middle panel, 
but do not include any upper limits in the right-hand panel. Com- 
pared to the whole sample (i.e., left-hand panel), there is a distinct 
lack of L1.4 GHz< 10^^ W Hz^i AGNs with very broad FWHMAvg 
and the rise in upper envelope with L[o m] is much less pronounced. 
This difference is highly unlikely to be due to the smaller numbers 
in the Li.4ghz< 10^"^ W Hz^' sample compared the whole sam- 
ple, with a KS-test indicating a ^ 0. 1 per cent chance that the two 
samples are drawn at random from the same underlying population. 

The right-hand panel of fig. [TO] shows the FWHMAvg versus 
Ljoiu] distribution for Li 4ghz> 10^"^ W Hz^'0 Despite there 
being far fewer sources in this plot overall compared to the mid- 
dle panel, over seven times more of them lie at extreme values of 
FWHMAvg (i-e., > 1000 km s"^; 59 compared to 8) and the en- 
tire distribution is centred around a higher FWHMAvg- We find that 
FWHMAvg shows no clear dependence on Lp in] for these radio 
luminous AGNs. 

In the left-hand panel of fig. [TT] we plot FWHMAvg against 
il.4GHz to further explore the connection between these two pa- 
rameters. We have included upper limits on Li 4 ghz which dom- 
inate the numbers at Li 4 ghz< 10^"* W Hz^'. The most striking 
aspect of this plot is the peak in the FWHMAvg distribution around 
^1.4 GHz= 10^"' — 10^^ W Hz^', which are the L] 4 ghz associated 
with the broadest of our stacked [O IIl]A.5007 profiles. As well as 
being traced by the detected sources, this peak is also defined by 
the upper limits on the radio-undetected sources. Furthermore, a 
KS test indicates that this distribution remains statistically the same 
if we restrict the sample to z < 0.123, at which point the flux limit 
of the corresponds to a maximum Li 4 ghz 

of 10^^ W Hz"', m- 

dicating that the distribution is not produced by a selection effect 
caused by the redshifts of the sources. At first sight this distribution 
suggests that moderate luminosity radio AGNs are most affective at 
producing broad [O IIl]A.5007 emission lines. However, other fac- 
tors must first be ruled out before we can arrive at this conclusion as 
it is possible that a combination of selection effects could produce 
this distribution. 

A possible explanation for the shape of the distribution at 
il.4GHz< 10^"' W Hz^' is that the correlation between Li,4ghz 
and L[o m] would make the broad component in these less radio 
luminous AGNs more difficult to detect by our line fitting proce- 
dure due to the weakness of the Ljq m] line. However, when we 
reproduce the FWHMAvg versus Li 4 ghz plot using AGNs with 
LjO iil]> 10"*^ srgs s^' to ensure that any broad component is de- 
tected (fig. [To] shows that any such detection bias is negligible at 
these [O IIl]A,5007 luminosities), the same distribution is produced. 
As such, we can rule out this detection bias as having a major in- 
fluence on the form of this distribution. 

The small numbers of AGNs in our sample with Li.4ghz> 
10^^ W Hz^' may be the reason for the apparent drop-off in peak 
FWHMavi. at these high Li 4 ghz- However, a KS test reveals that 
this is unlikely to be the case, returning only a --^0.05 per cent prob- 
ability that the FWHMAvg of L1.4 ghz= 10^"^ - 10^^ W Hz"' and 
are drawn at random from the same underlying 
distribution. This and the arguments laid out above lead us to con- 
clude that the peak in the FWHMAvg-^'i.4 GHz distribution around 
^I.4GHz= 10^"^ ~ 10^^ W Hz^' is unlikely to be due to selection 

* We note that the numbers of AGNs in the middle and right-hand panels 
do not sum to the numbers in the left-hand panel as there are significant 
numbers of AGNs in our sample with Li 4 ghz upper limits > 10^^ W Hz"' 
which could not be reliably included in either the middle or right-hand plots. 
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Figure 10. The distribution of FWHMAvg (see Eqn.|2) plotted as a function of L|q mi for all (left-hand panel), Li 4 gHz< 10^' W Hz ' (middle panel) and 
^■1.4 GHz> 10"' W Hz"' (right-hand panel) AGNs. The upper envelope of this distribution peaks around lO"*^ ergs s"'. However, the comparing the middle 
and right hand panels reveals that the AGNs with the broadest tend to be the more radio luminous of our sample. 



effects and thus gives a true impression of the relationship between 
the profile of the [O III] A.5007 emission line and the radio properties 
of the AGN. 



5 DISCUSSION 

In the previous sections we have explored how [O III] A.5007 emis- 
sion line profiles differ between AGNs of different Lp^Qfiifxom. 
^[Olll|)' ^1.4 GHz! ^Edd™d '^Rad- Since emission line profiles are 
governed by the kinematics of the emitting gas, our results allow us 
to assess how the movement of gas surrounding the AGNs is related 
to these four key observables. In this section, we place our results in 
the context of previous studies of the kinematics of the [O IIl]A.5007 
emitting gas and, in doing so, consider the implications of our re- 
sults on our understanding of how AGN-driven outflows influence 
galaxy growth. 

The results from both our spectral stacks and our individual 
[O IIl]A,5007 line fits shows that, of the four observables consid- 
ered (i.e., Lagn> 1-\A GHz. ^Eddand Rr^a), Lj 4 ghz is the most im- 
portant factor in dictating the profile of this emission line. Our 
results indicate that the width of the [O III] ^5007 line peaks be- 
tween Li.4GHz= 10^^ W Hz-i and ~ 10^^ W Hz"' with both 
more and less radio powerful AGNs having, on av erage, narrower 

t o III A.50 07 lin es. Similar results were reported bv lHeckman et al.l 
19841) and|whittle ( 1992) for much smaller samples of AGNs. The 
size and relative completeness of our sample allows us to build on 
these previous studies, breaking the degeneracy between Ljq ni| 
and ii. 4 GHz to identify the underlying links, and placing AGNs 
with kinematically energetic narrow line regions in the context 
of the general population by quantifying the fraction with broad 
[O IIl]A,5007 emission lines. 

Spectral studies of small samples of AGNs displaying com- 
pact radio jets have reported broadened [ O IIl]A,5007 emission 
lines that are co-spatial with these jets (e.g.. lTadhunter et"al]|2003l : 
iHolt et aL 2006, 2008), suggesting that it is the jets' mechanical 
energy that is disturbing the [O IIl]A,5007 emitting gas. Since we 
have no information regarding the spatial extent of the [O III] A.5007 
emitting gas in our sample, we are unable to perform such a 
direct analysis here. Instead, we can exploit the spatial resolu- 



tion of the FIRST radio survey to determine whether AGNs with 
broad [O IIl]A,5007 have predominantly extended or compact radio 
morphologies. Taking NVSS-detected AGNs in our sample with 
il.4 GHz= 3 X 10^3 - 3 X 10^4 W Hz-1 (to cover the peak of the 
FWHMAvg distribution), we extract the angular extent of their ra- 
dio major axis from the cross-matched FIRST catalogue, which 
has a higher spatial resolution at 1.4 GHz compared to NVSS. We 
find all AGNs in this sample with FWHMAvg> 1500 km s"' and 
single FIRST matches have deconvolved extents < 2" (i.e., unre- 
solved) and there is no evidence that more extended radio sources 
have broader [O IIl]A,5007 lines. Furthermore, in the case of those 
with multiple FIRST matches, the vast majority of AGNs in this 
subsample (i.e., > 80%) have at least 1 radio component closer 
than 1" to the galaxy nucleus, suggestive of a radio core. Finally, 
none of the AGNs in this radio-selected subsample have multi- 
ple NVSS matches, meaning they are unlikely to have extended, 
lobe-dominated radio emission without also having a compact core. 
While not conclusive, these various outcomes support the view that 
compact (i.e., ~kpc-scales or less) radio cores play a major, if not 
the dominant, role in strongly disturbing the [O IIl]A,5007 emitting 
gas around optically-selected AGNs. This may also explain why 
the [O III] A.5007 width peaks at moderate Li 4 ghz. close to the di- 
vision between Faranoff-Riley types I and II, since the latter dissi- 
pate more of their energy at larger (i.e., 3> kpc) scales. However, 
we also point out that the peak of the FWHMAvg distribution is 
close to the point at which the radio luminosity function becomes 
AGN, rather than star-formation, dominated. It is therefore feasi- 
ble that the radio output of some of the AGNs with the broadest 
[O IIl]A,5007 lines is dominated by star-formation. 

One factor that we have not yet explored is the possibility 
that some of the AGN in our sample may have boosted radio 
emission due to beaming, meaning that the intrinsic radio power 
may be considerable lower than what we infer here. However, 
while this may explain part of the decreased disturbance toward 
higher Li 4 ghz> it is unlikely that all (or even the majority) of the 
^1.4 GHz^ 10^^ W Hz"' AGNs in our sample are beamed (i.e., 
blazars), especially since emission line selection would tend to bias 
against such AGNs as they tend to be continuum dominated (e.g., 
iMorris et al ]|l99lUStockeetal.ll99ll) . 

Our analyses show a higher fraction of moderate radio- 
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Figure 11. Distribution of FWHMAvg plotted as a function of Li 4 ghz- The left-liand panel includes all of the AGNs in our sample, whereas the right-hand 
panel show only those with Lp iii|> 10**^ ergs s^' . Both panels show a significant (i.e., 0.05 per cent chance of being produced by random) peak in the upper 
envelope of the distribution around Li 4 gHz~ lO^"* W Hz"' . 



luminosity (i.e., Li.4ghz= lO^"' - 10^^ W Hz"') AGNs display 
broad [O IIl];^5007 lines compared to either weaker or stronger 
radio AGNs. Furthermore, these fractions are only weakly depen- 
dent on Lagn- These moderate Li 4 ghz AGNs represent roughly 
3 per cent of optically-selected AGNs at z < 0.123 (the redshift at 
which the NVSS survey is complete to Li 4 ghz> 10^^ W Hz"'). 
As such, while they constitute a small fraction of optically-selected 
AGNs overall, they are at least ~10 times more common than the 
more radio luminous AGNs that are typically the focus of detailed 
studies of AGN outflows (e.g.. iTad hunter et al. 2003; Holt et al. 
20031: iNesvadba et aL.2003 ; lHolt et al. 2008.^Nesva dba et al. 2008. ; 



Cano-Dfaz et alj2012l ; see Fig. 1 of lHarrison et al.l20l'2l fora com 

parison). Those studies that have sampled the Li,4ghz~ 10^'' — 
10^^ W Hz"', i[om]> lO'*^ srgs s"' range have reported out- 
flows extended over galaxy (i.e., ~few kpc) sc ales, suggesting that 
they could be affectin g their host galaxies (e.g.. lGreene et al .I2OIII ; 
iHarrison et al.ll2012l) . A detailed comparison of the resolved ra- 
dio components and the gas kinematics is required to determine 
whether (and, if so, how) the two are physically linked. 



6 SUMMARY 

We have used a large (i.e., 24264) sample of optically selected 
AGNs from the SDSS database to explore how the profiles of the 
forbidden [O IIl]A,5007 emission line relate to other key AGN pa- 
rameters (bolometric luminosity, Eddington ratio, radio loudness 
and radio luminosity) . We use spectral stacking analyses to deter- 
mine how the average [O IIl]X5007 profile changes as a function of 
these four parameters . We also explore the fractions of AGNs in 
each of these bins that have broad [O IIl]A,5007 emission lines. Our 
main results can be summarised as; 

• The average [O ni]A.5007 profile of Type 1 AGNs displays a 
prominent blue wing that can be well described by a broad (i.e., 
FWHM= 851 km s"') gaussian component that is blueshifted by 
148 km s"' from a narrow (i.e., FWHMl = 335 km s"') core. The 



flux contained within this broad component represents 45 per cent 
of the total average line flux (see i)3.1t . 

• When we separate the Type 1 AGNs into discrete bins of 
L[o III] or Eddington ratio (A-Edd) we find that the blue wings on 
the average [O ni]A.5007 profiles are increasingly prominent at 
higher values of Ljq m] ^nd in our highest A-Eddbin (i.e., A,Edd> 0.3; 
see i)3.2.1l and fig. |3]l. However, when we correct for the known 
correlation between i[o m] ^'1.4 GHz we find that Li.4ghz 
has the most profound effect on the [O IIl]A,5007 profile, with 
L\A GHz= 10^^ - 10-^ W Hz"' AGNs having the broadest profiles. 
This appears to be the case for high Eddington ratio AGNs (i.e., 
^Edd> O-l)' although we are not able to show this conclusively in 
our highest Xgdd bin (A,Edd> 0.3) due to the small number of such 
extreme AGN even in our large sample. 

• There is no difference in the average [O ni]A.5007 profiles of 
broad-line and narrow-line Seyfert 1 AGNs. 

• When the [O IIl]A,5007 profiles of Type 2 AGNs are stacked in 
terms of Ljg m] their [O IIl]A,5007 widths also increase with Ljq m] 
but are more symmetric than those of Type 1 AGNs. We interpret 
this as suggesting that the blue wings are the result of outflowing 
material directed roughly along the axis of the obscuring torus. 

• Overall, roughly 25 per cent of all AGNs have average 
[O III]A,5007 FWHM (i.e., FWHMAvg) > 500 km s"', compared 
to -50 per cent of AGNs with L1.4 ghz> 10^^ W Hz"' . All of the 
AGNs in our sample with FWHMAvgSi 1200 km s"' have radio 
detections or upper limits consistent with L1.4 ghz> 10^"' W Hz"' . 



Acknowledging the limiting spatial resolution of the FIRST 
and NVSS surveys, our results are consistent with models in which 
compact radio cores play an important role in disturbing the dif- 
fuse gas around the AGN. Ascertaining whether the impact of these 
cores is felt on extended, galactic scales will require spatially re- 
solved radio and spectroscopic observations. 
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APPENDIX A: MULTI-COMPONENT FITTING ROUTINE 

The emission line fitting procedure we have developed starts with a 
simple, single gaussian fit to the [O IIl]XX4959, 5007 emission lines 
and progresses by increasing the complexity of the fit until each line 
is fit with multiple gaussians where necessary. At each stage of the 
fitting procedure any change in x^is evaluated to establish whether 
the addition of any new gaussian components has resulted in an 
improved fit (at a >99 per cent confidence level). 

First, the procedure attempts to fit each of the 
[O IIl]XM959,5007 emission lines simultaneously with sin- 
gle gaussians. We require the FWHM and velocity offsets of 
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these gaussians to be equal throughout the fitting process and the 
normahsation of the [O IIl]A.4959 gaussian fixed to 1/3 that of the 
[O III]?l5007 gaussian (i.e., FWHMp ni];v4959 = FWHMp m]A,5007, 

l'[0 III]?.4959 = ^[O III]A,5007' -^[O III]?t4959 = 1/3 ' ^[O III]?t5007)- By 
fitting [O IIl]A.4959 we minimise any influence both this and any 
neighbouring emission lines (e.g. the Fe II complexes) may have 
on the measurement of the [O IIl]A,5007 profile. We note, however, 
that we do not use the results from fitting the [O IIl]X4959 in our 
analysis. The routine then attempts to fit each of the permitted HP 
and Ha lines simultaneously with single gaussians (referred to later 
as the 'original' HP and Ha gaussians), the velocity offsets and 
FWHM of which are, again, linked (i.e., FWHMep = FWHMhk, 
Vjjjj = v'Ha)- Second gaussians are then added to the fits of the 
[O III] A,X4959, 5007 lines in an attempt to model any complexities 
in the line profiles. The velocity, FWHM and normalisation of 
this second set of [O III] A,A.4959, 5007 gaussians are linked in the 
same manner as the first set. Finally, the procedure attempts to 
fit the narrow components of the HP and Ha lines and each of 
the [N Il]A,A,6548,6584 lines with up to two further gaussians, de- 
pending on whether the addition of the second [O iii]XX4959, 5007 
gaussians resulted in a significantly improved fit. Each of these 
additional gaussians are forced to have the same velocity offsets 
and FWHM as their respective [O IIl]A.A.4959,5007 counteiparts, 
but their relative normalisations are allowed to vary. We caution, 
however, that it has previously been reported that t he line profiles 
of d ifferent forbidden species can differ (e.g., iKomossa et al.l 
12008) . By adopting the approach that the relative normalisation 
of the two [O IIl]A.5007-derived narrow-line components can 
vary, we adopt a compromise between having some reasonable 
constraints on the fit to the permitted and [N Il]X6584 lines (thus 
preventing these lines from broadening to accommodate the true 
broad permitted lines), while allowing the profile to vary somewhat 
to accommodate some of the differences between different species. 
The normalisation of the [N Il]^X6548 line is fixed to 1/3 that of 
the [N Il]A.6584 line. In fitting the narrow components of the HP 
and Ha lines with the modelled [O III] XX4959, 5007 profile, the 
original HP and Ha gaussians are allowed to widen to fit any broad 
components of these permitted lines. 

Since the goal of our fitting the emission line profiles is to 
look for general trends within the large AGN population, we are not 
greatly concerned if a small number of "one off" spectra are poorly 
fit by our routine, provided that it provides accurate fit parameters 
for the vast majority of sources in our sample. To demonstrate that 
this is the case, we plot the results of stacking randomly selected 
spectra after subtracting the best fit model (see the red line in fig. 
lAlb . We perform this for separate samples of BLSls, NLSls and 
Type 2 AGNs, randomly selecting 1000 AGN from each of these 
three samples. In general, we see only weak residuals around the 
fitted emission lines, especially the [O IIl]X5007 lines (the bump 
around the [O III] A.5007 line in the NLS 1 average spectrum is due to 
the Fe II complex which is stronger, in general, in this class AGNs). 
There is evidence of a narrow residual in the Ha and HP lines of 
the NLSls (with a peak roughly 1/lOth the maximum flux of these 
lines). However, since we only use the broad components of these 
lines to discriminate between NLSls and BLSls, which are gener- 
ally well fit, this will not affect our main results. As such, we take 
these average residuals as evidence that, in general, our line fitting 
routine is providing an accurate and unbiased parameterisation of 
the emission profiles of our sample spectra. 
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Figure Al. Stacked emission line profiles (left: Hp, [O IIl]^5007; right: Ha) for 1000 BLSls (top), NLSls (middle) and Type 2s (bottom) AGNs randomly 
selected from our full sample. The red lines show the results of stacking the spectra after subtracting the best fit line model (i.e., — We 

take the general lack of strong systematics in these average residual spectra as evidence that, over the sample, our line fitting routine is doing a reasonable job 
of modelling the profiles of the emission lines, particularly the broad wings of the Hct and Hp fines and the overall [O IIl]^5007 profile. 



